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Abstract 



We have studied the null geodesies of the Schwarzschild black hole surrounded by 
quintessence matter. Quintessence matter is a candidate for dark energy. Here, we 
have done a detailed analysis of the geodesies and exact solutions are presented in 
terms of Jacobi-elliptic integrals for all possible energy and angular momentum of the 
photons. The circular orbits of the photons are studied in detail. As an application 
of the null geodesies, the angle of deflection of the photons are computed. 
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1 Introduction 

It is well know that we live in a universe undergoing a period of accelerating expansion. 
This indicates the presence of mysterious form of repulsive gravity called dark energy. 
There are several observations that support dark energy. One is the studies of type 
la supernova [1] [2]. The other is the observations related to the cosmic microwave 
background (CMB) [3] and large scale structure (LSS) [4] [5]. 

The nature of dark energy is yet to be understood. There are several cosmological 
models proposed in which the dominant component of the energy density has negative 
pressure. One of them is the cosmological constant A which corresponds to the case 
of dark energy with a state parameter Uq = —1. However, there is a key problem 
that is yet to be understood about the cosmological constant from the point of view 
of fundamental physics. Its observed value is too small which is also called the fine- 
tuning problem [7]. 

There are alternative models that are proposed as candidates for dark energy. 
Most of these models are based on a scalar field. Such scalar field models include 
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but not limited to, quintessence [8], chameleon fields [9], K-essence [10], tachyon field 
[11], phantom dark energy [12] and dilaton dark energy [13]. For a detailed review of 
various models of dark energy see [6] . 

Given the fact that dark energy contains about 70% of the universe and black holes 
are also accepted as part of our universe, studies of black holes surrounded by dark 
energy takes an important place in research. In this paper, we study the Schwarzschild 
black hole surrounded by quintessence matter. Quintessence is described as a scalar 
field coupled to gravity with a potential that decreases as the field increases [8] [6]. 
Kiselev [14] derived black hole solutions surrounded by quintessence matter which 
has the state parameter in the range, — 1 < 0;^ < — |. In this paper, we will focus on 
the solution derived by Kiselev. 

There are several works done on the black hole derived by Kiselev in the literature. 
The quasinormal modes of the black holes has been computed extensively. The quasi- 
normal modes of the Schwarzschild black hole surrounded by the quintessence matter 
has been computed in [15] [16] [17] [18]. The quasinormal modes of the Reissner- 
Nordstrom black hole surrounded by the qiiintcssence matter has been presented in 
[19]. Hawking radiation of d-dimensional extension of the Kiselev black hole has been 
studied by Chen et. al. in [20]. 

The main objective of this paper is to study the geodesic structure of massless 
particles of the Schwarzschild black hole surrounded by quintessence matter. The 
motion of particles around black holes has been studied extensively in the literature 
for all types of black holes. Due to the large volume of papers related to this subject, 
we will avoid referring the papers here. Motion of p [articles around a black hole is one 
way to understand the gravitational field around a black hole. Given the fact that 
dark energy is one of the most important issues need to be resolved in physics, the 
studies of particles around a black hole immersed in dark energy take an important 
place. 

The paper is presented as follows: In section 2 the black hole solutions surrounded 
by quintessence matter is introduced. In section 3 the general formalism of the 
geodesies are given. In section 4, the radial geodesies are discussed. In section 5, 
the geodesies with angular momentum are discussed. In section 6, the geodesies are 
analyzed with a new parameter. In section 7, the applications of the null geodesies 
are given. Finally, the conclusion is given in section 8. 

2 Schwarzschild black hole surrounded with 
quintessential matter 

In this section we will give an introduction to the Schwarzschild black hole solutions 
surrounded by quintessential matter obtained by Kiselev[14]. The geometry of such 
a black hole has the metric of the form, 

ds"" = -g{r)de + g{r)-^dr'' + r^dO'' + sin''ed(t>^) (1) 
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where, 

2M c 

9^')-^- — -^:^. (2) 

Here, M is the mass of the black hole and Wq is the state parameter of the quintessence 
mater, c is a normalization factor. The parameter Wq has the range, 

-l<wq< -i (3) 
The equation of state for the quintessence mater is given by, 

Pq = WqPq (4) 

and 

_ c 3Wq 

Pi- 2 r^{i+w,) 

Here Pq is pressure and Pq is the energy density. As described in [14], to cause 
acceleration, the pressure of the quintessence matter Pg has to be neagtive. The 
mater energy density Pq is positive. Hence the parameter c > for negative Wq. For 
more details on the derivation of the solutions and the basis to choose the parameters 
as given, reader is referred to the paper by Kiselev [14] . 

One can observe that for ty^ = — 1, the function g{r) for the metric reduces to 

, . 2M 2 
g{r) — 1 cr (6) 



which is the Schwarzschild-de-Sitter black hole space-time 

^9 ~ ^3 



In this paper, we will focus on the special case of Wq — Then, 



/ N 2M 
9{r) = 1 cr (7) 

For 8Mc < 1, the metric with the above g{r) given in eq.(7) has two horizons as. 



1 - VI - 8Mc 
rin = (8) 

1 + VI - 8Mc 
rout = ^ (9) 

The inner horizon is like the Schwarzschild black hole horizon. The outer horizon 
Tout is a cosmological horizon similar to what is observed in the Schwarzschild-de- 
Sitter black hole. Notice for small c, r„^t ~ ^ which is similar to Tout ~ for 
small A in the Schwarzschild-de-Sitter case. There is a static region between the two 
horizons. Hence, in considering the motion of photons, we will focus on the region 
between the two horizons. 

If SMc — 1, the space-time has a degenerate solution at r = ^. For 8Mc > 1, 
there are no horizons. All three situations are represented in Fig. 1. 
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Figure 1. The graph shows the relation of g{r) with r for various values of the pa- 
rameter c. Here, M = 1 

The Hawking temperature for this black hole is, 



An 



dgtt 



dr 



(10) 



which leads to, 



T 

■J- 7 



1 

An 



2M 



^in,out 



(11) 



The temperature for both horizons are plotted in Fig. 2. The temperature of the outer 
horizon is smaller than the inner horizon. 



Temp 0.005 - 
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Figure 2. The graph shows the relation of the temperature with c for both horizons. 
The dashed curve is for the inner horizon while the thick curve is for the outer hori- 
zon. Here, M = 2 
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3 Geodesies of the black hole 




We will derive the geodesic equations for neutral particles around the quintessence 
black hole. We will follow the same approach given in the well known book by 
Chandrasekhar[21]. The equations governing the geodesies in the space-time given in 
section(2) can be derived from the Lagrangian equation, 

^ - - 1 -sM + ^ g)' + (f )^ ^ ISM (^2) 

Here, r is an affinc parameter along the geodesies. Since the black hole have two 
Killing vectors dt and d^p, there are two constants of motion which can be labeled as 
E and L. They are given by, 

g{r)i = E (13) 

r^sin^e^) = L (14) 

We choose 9 — and ^ = as the initial conditions. Hence 9 — 0. 9 will remain 
at 7r/2 and the geodesies will be described in an invariant plane &t 9 — 7r/2. With i 
and 4> given by equations(13) and (14), the Lagrangian in eq.(12) simphfies to be. 

Here, 2£ = h. h = 1 corresponds to time-like geodesies and h = corresponds to 
null geodesies. For a time-like geodesic, r may be identified with proper time of the 
particle describing the geodesic. Comparing eq.(15) with -|- Veff — E^, one get the 
effective potential, 

Veff =[^ + h) 9{r) (16) 

By eliminating the parameter r from the equations(13) and (14), one can get a relation 
between and r as follows; 

# L 1 ^^^^ 



In this paper, we will only focus on null geodesies. Hence h = and the effective 
potential corresponds to. 



Veff = L^'-^ (18) 
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4 Radial null geodesies 



The radial geodesies corresponds to the motion of particles with zero angular mo- 
mentum (L = 0). Hence for radial null geodesies, the effective potential, 

Veff = (19) 

The two equations for i and r simplifies to. 



f = ±E; 

The above two equations lead to, 

dt 1 
dr g{r) 



E 
9{r) 



2M 



(20) 



(21) 



(1 - =^ - cr) 

The above equation can be integrated to give the coordinate time t as a function of 

r as, 

2 „ , _i / -1 + 2cr ^ 



t = ± 



zTanh 



Log (r - 2M 



cr 



+ const± (22) 



Vl-8cM VVl -8Mc, 

The constant of integration is imaginary. However, when one compute the time t, the 
overall expression will be real. On the other hand, the proper time can be obtained 
by integrating. 



dr 
dr 



±- 



'E 



which leads to, 



r = ±— + const ± 
E 



(23) 



(24) 



When r ^ rin, t 



oo. When r — )■ r^^, r — )■ which is finite. Hence the proper 
time is finite while the coordinate time is infinite. This is the same result one would 
obtain for the Schwarzschild black hole[21]. Fig.3 represents both times to show this. 
Note that we have picked the plus sign in this case, since we are studying the ingoing 
light rays. 
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Figure 3. The graph shows the relation of the coordinate time t and the proper time 
r with the radius r. Here, M = 1, c = 0.05. The initial position is at r = 16. The 
inner horizon is at rj„ = 2.254 and the outer horizon Tout = 17.746. The dark curve 
is for t and the light curev is for r. 

5 Null geodesies with angular momentum (L ^ 0) 

In this section, we will study the null geodesies with angular momentum. 

5.1 Effective potential 

In this case, 

Veff = 9{r)^ (25) 

For r = Tin and Tout, Veff = 0. In Fig. 4, the Veff is given for various values of 
c. The potential for the Schwarzschild black hole is higher than for the black hole 
surrounded with the quintessence. 
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^^^^^^^^^^^^^^^^^ c = 0.05 







Figure 4. The graph shows the relation oi Veff with the parameter c. Here, M = 1 
and L = 1 

When the Veff is expanded, it looks like, 

^^fS = ^-^^ (26) 

The contribution to the effective potential from the quintessence is — — compared 
to the Schwarzschild black hole. Note that in the Shwarzschild-de-Sitter case, the 
contribution is a constant — L^A. 
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Since the energy and the effective potential are related by + K// = E^, the 
motion of the particles depend on the energy levels. In Fig. 5, the effective potential 
is plotted and three energy levels, Ei, Ec and E2 are given. There are three different 
scenarios of motion of the particles as described below. 

Case 1: E = E^ 

Here, E"^ — Veff = and r = 0. The orbits are circular. Due to the nature of the 
potential at r = r^, these are unstable. More details of the circular orbits will follow 
in the section(5.2). 

Case 2: E = E2 

Here, E^ — Veff > in two regions as is clear from Fig. 5. If the photons start the 
motion at r > r^, it will fall to a minimum radius and fly back to large values of r. If 
the photons start the motion at r < r^, then the photons will fall into the singularity 
crossing the inner horizon at r = rj„. 

Case 3: E = Ei 

Here, E"^ — Veff > and r > for all r values. Hence the photons coming from 
large r values will cross the horizon at r = rj„ and will fall into the singularity. 
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Figure 5. The graph shows the relation of Veff with the energy E. Here, M = 1, c = 
0.01 and L = 1. 

5.2 Circular orbits 

Circular orbits occure when E 

and 
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= Ec as explained in section(5.1). Hence at r = rc, 
r = O^Veff = El (27) 



^ = (28) 
dr 



Prom eq.(28), one get two solutions for circular orbit radius r as, 



r± = ^ (l ± Vl - 6cM) (29) 

The larger root r+ is greater than rout and r^n < r_ < rout- Since the motion of 
interest is between the horizons, the radius of circular orbits considered will be r_. 
We will name it as Tc in the rest of the paper. Hence, 

= - (l - Vl - 6cM) (30) 



' c 

C 



Due to the nature of the potential at r = Tc, the circular orbits at r = Tc are unstable. 
The hypersurface at r = is known as the "photon sphere" . For a detailed discussion 
about photon spheres see the paper by Claudel et.al[22]. When c — )■ 0, Tc — )■ 3M 
which is the radius of the unstable circular orbit of the Schwarzschild black hole[21]. 
Circular orbits takes a special place in the studies of geodesies. The circular orbit 
is the boundary between two qualitatively different regions. An interesting paper on 
null circular orbits of a hairy black hole is given by Hod [23]. 

The radius of the circular orbit in eq.(30) is independent of E and L. However, 
they are related to each other from eq.(27) as, 

E', g{T.) _ {r.-2M-crl) 1 

LI rl rl Dl ^ ^ 

Here, Dc is the impact parameter at the critical stage. When c — > 0, — > 27M^ 
which is the impact parameter for the unstable circular orbits of the Schwarzschild 
black hole [21]. 

5.2.1 The Time Period 

The time period for the circular orbits for proper time (r) as well as coordinate time(i) 
can be calculated with = 27r. From eq.(14), 

27rr^ 

t; = (32) 

Prom combining the two equations eq.(13) and eq.(31), 

„ 27rrc 27rry2 



rc-2M- crl 



(33) 



The values and Tt for the Schwarzschild black hole is given by, 

Tt,sh = 3\/3M (34) 
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and 

T... ^ (35) 

The graphs for the time periods are given in Fig. 6 and Fig. 7. By observing the graphs 
of the time periods, it is clear that the periods for the Schwarzschild black hole is 
smaller in comparison with the quintessence black hole. In a recent paper, Hod made 
an interesting observation that the null circular geodesies provide the shortest possible 
orbital period to circle a black hole [24] . 



T(piopei" time) 




M 



Figure 6. The graph shows the relation of T^. with the mass M. The dark curve is 
for Tt- ^quintessence and the dashcd curve is for T^-^sch- Here, c = 0.05 and L = 1. 



T(cooridnate time) 




M 



Figure 7. The graph shows the relation of Tj with M. The dark curve is for 
Tt,quintessence and the dashcd curve is for Tt^sch- Here c = 0.01 and L = 1. 

5.2.2 Lyapunov exponent for the unstable circular orbits 

The instability time-scale of the unstable circular null geodesies is given by the Lya- 
punov exponent A. The expression for A was derived by Cardoso et.al [25] as. 
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Figure 8. The graph shows the Lyapunov exponent A as a function of c. Here, M = 1. 

The instabihty of the circular orbits are more for the Schwarzschild black hole in 
comparison with the black hole with a non-zero c. 

5.3 Force on the Photons 

Since we have already computed the effective potential, one can obtain the effective 
force on the photon as, 

2 dr r3 2r^ ^ ' 

Here, we have divided by 2 since the equation of motion is written as + 

Veff = E"^. The first term, which is the Newtonian term, is attractive since it is 
negative. The second term is repulsive. The third term which is the force due to the 
quintessence matter, is attractive. It is interesting to notice that even though dark 
energy in cosmology is associated with a repulsive force to facilitate acceleration, in 
this situation, the dark energy term is attractive. This is a result due to the fact that 
we are studying a static configuration between the two horizons. 
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Figure 9. The graphs shows the relation of F with r. Here M = 0.5, c = 0.1 and 
L = 2. 

Prom Fig. 9, the force F = at, 

r± = - (l ± Vl - 6cM) (38) 



One can see from the graph that rout < (l + Vl — 6cM) . Since we are only concern 
about the motion for r < Tout, one can conclude that the force is positive for < 
< ^out which leads to a repulsive force. On the other hand, for r < r^, the force is 
negative leading to an attractive force. 

The maximum repulsive force occur at r = r^ax where, 

6 Analysis of the geodesies with the variable u = ^ 

Geodesies equation of motion can also be studied using a well known change of variable 
u = -. Then the eq.(17) becomes. 
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where, 



f{u) = 2Mu^ - + CU + 



1? 



(41) 



It is clear that the geometry of the geodesies depends on the nature of the roots 
of the equation f{u) = 0. Note that for any value of the parameters in the theory 
M,c,L,E, the function f{u) — )■ — oo for u — )■ — oo and f{u) — )• +oo for u — > +oo. 
Also, for u = 0, f{u) = +j2- Therefore, f{u) definitely has one negative real root 
{ui). Let us name the roots of f{u) as Mi,U2 and M3. The sum and the products of 
the roots mi,M2 and M3 of the polynomial f{u) are related to the coefficients of f{u) 
as [26], 

Ml + M2 + Ms = ^ (42) 

"^"^"^ = -2ML^ ^^^^ 
As discussed earlier, ui is real and negative. Therefore, the roots U2,u^ has to be 
positive if they are real. This conclusion comes from observing the signs of the eq.(42) 
and eq.(43). Overall, the polynomial f{u) has a negative(tii) real root always. The 
other two roots (^2, W3) will be either real or complex-conjugate. If they are real, then 
they will be positive. Also, if they are real, they could be degenerate roots as well. 
The function f{u) for general values of M, c, E and L is given in Fig. 10. 




Figure 10. The graph shows the function f[u) for M 
E = 9. 



0.5, c = 0.01, L = 30 and 



6.1 General solution for the geodesies in terms of u 

From eq.(40) derived above. 
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where f{u) can be written as, 

f{u) — 2M{u — ui){u — U2){u — u^) 



(44) 



(45) 



The "+" sign will be chosen without lose of generality. One can integrate the equation, 
—p= = dcf) to get a relation between u and in terms of Jacobi-eUiptic integral 

= J ^^'^-^ _|_ constant (46) 



yj2m{u2 - Ui) 



Here, 



sin^ 



y = 



{U2 - Ui) 

(u-ui) 

{us - Ui) 



(47) 



(48) 



{U2 - Ui) 

Depending on the values of the E, L, c, M, the integration constant would be real or 
imaginary. Also, Ui is always real while the nature of U2 and depend on the values 
of E, L, c, M. 



6.2 Circular orbits 

We will analyze in more ddetail the circular orbits here. In this case, the function 
f{u) has a degenerate root at u — Uc- Hence, f{u) can be written as, 

f{u)^2M{u-u^f(u-ui) (49) 

Here, Uc — ^i where Tc is given by eq.(30). From eq.(42), 

1 / ^ 

u^^ - 2uc 50) 

A null geodesic arriving from r — rs > Tc { or u — Us < Uc) undergo unstable circular 
orbits at r — rc { or u — Uc)- The from of the function f{u) in this case is given in 
Fig.ll. 
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Figure 11. The function f{u) when it has degenerate roots leading to circular or- 
bits. The shaded area represents the possible regions of motion. Here, M = 0.5, c = 
0.01 and r, = 1.51142. 



One can integrate the equation, p— = dcf) to get a relation between u and as. 



u = Ui + (uc — Ui)tanh (51) 



bo and are given by. 



bo = aotanh Ma/— -) (52) 

\V Ui-UcI 



V M{Uc - Ml) 

In Fig. 12, the polar plot of the null geodesies are given for photons arriving from 
Ts > Tc and having an unstable circular orbit at r = Tc. 
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Figure 12. The polar plot shows the critical null geodesies approaching the black 
hole from > r^. The geodesies have an unstable circular orbit at r = r^. Here, 
M = 0.5, c = 0.01 and = 1.51142. 

6.3 Unbounded orbits 

Case 2 : E = E2 

Here, we will study Case 2 given in the section(5.1). In this case, f{u) = has 
three real roots. Therefore, the motion is possible in two regions given as 1 and 2 in 
Fig 13. 



f(u) 
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Figure 13. The graph shows the function f{u) when it has three real roots. Here 
M = 0.5, c = 0.01, L = 30 and E = 9 

If the particle starts far from the black hole at r = r^, it will fall until r = r2 = l/u2 
and fly away from the black hole. The equation of motion is the same as given 
in eq.(46). The integration constant is chosen such that = for r = r^. The 
corresponding motion is given in Fig 14. 
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Figure 14. The polar plot shows the null geodesies approaching the black hole from 
r = r,. Here, M = l,c = 0.01, L = 30, ^ = 9,r, = 15 and = 2.54082. The inner 
horizon is at rj„ = 1.01021 

The second possibility for E = E2 corresponds to the motion starting from r = 
(or M = M3 ). Here, < Vc- Hence the particle will fall into the singularity crossing 
the horizon at r = rj„. In this case, the solutions for (p is similar as in eq.(46). 
Therefore, we will omit the explicit expressions for the in this case. The integration 
constant is chosen such that = when u = u^. The corresponding motion is given 
in Fig 15. 
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Figure 15. The polar plot shows the null geodesies falling into the black hole from 
r = rg. Here, M = 1, c = 0.01, L = 30,E = 11, rg = 1.34587 and = 1.01021. 

Case 3 : E = El 

Here, we will study the Case 3 given in the section(5.1). In this case, f{u) = 
has only one real root as given in Fig 16. 



f(u) 







/ 
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Figure 16. The graph shows the function f{u) when it has only one real root. 
Here, M = 0.5, c = 0.01, L = 21,E = 9 and Ui = -0.363391. 
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The motion is possible in the shaded area. For all values of r, the photons will 
fall into the black hole. The equation of motion is same as in eq.(46). However, in 
this case, U2 and us are both imaginary and ui is real. The corresponding motion is 
given in Fig 17. 



1.5 - 



1.0 - 



0.5 - 



0.0 



-0.5 - 



-1.0 - 



-1.5 - 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Figure 17. The polar plot shows the null geodesies approaching the black hole from 
r = 15 . Here, M = 0.5, c = 0.01, L = 21 and E = 9. 

7 Application of null geodesies: Gravitational lens- 
ing by the quintessence black hole 

Bending of light and gravitational lensing is an important aspect one can study re- 
lated to null geodesies around black holes. The unbounded orbit with energy E = E2 
studied in section(6.3) clearly represents how photons deviate from its original path 
when moving around a black hole. The motion is given in Fig. 14. 

7.1 Closest approach r^, 

To study the bending of light, let us first calculate the closest approach distance Vo 
for the photon with energy E2. It is defined by the value of r when ^ = 0. From 
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eq.(40) and eq.(41), 



/ 1 drV ^, , 2M 1 c E'' 

[-.d^) -nr)-^--. + - + J^ (54) 
^ = corresponds to the roots of f{r) = which also could be written as, 

+ D^cr^ - D\ + 2MD^ = (55) 

From Fig.l3, the roots of the above equation corresponds to ri,r2,r^ which are the 
inverse of Ui, U2, u^. Since ri < 0, the roots to be considered are r2 or ra. From Fig. 5, 
r2 and are the positions where E2 — Vgff. Clearly, < r2 (since U3 > U2)- Also, 
< Tc < r2- Hence, the root we will choose for the closest approach is r2- 
Now, one can use the well known techniques in determining the roots of a cubic 
polynomial to obtain tq from eq.(55) as. 



-■2— = 2ff-H^"'^"'(if?))-^ (56) 

Here, p and q are given by, 

^_P^c^ (57) 

54:MD'^ + 9cD^ + 2D^c^ , , 

? = (58) 

When c — )■ 0, To approaches the well known value for the Schwarzschild black hole 
given by [27], 

Fig. 18 shows the closest approach for the quintessence black hole and the Schwarzschild 
black hole. 
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Figure 18. The graph shows the closest approach for the quintessence black 
hole(dark curve) and the Schwarzschild black hole(dashed curve) as a function of the 
impact parameter D. Here, M = 1 and c = 0.05. Dc = 6.8956 for the quintessence 
black hole and Dr = 5.19615 for the Schwarzschild black hole. 



7.2 Bending angle 

Now, we will compute the angle of deflection of light for the quintessence black hole. 
In a paper by Amore et.al.[28], a new method was introduced to compute the deflec- 
tion angle for a static spherically symmetric space-time. This method yields highly 
accurate analytical results. In the paper, they applied the technique to well known 
metrics. One of them was the metric in Weyl gravity [29] [30]. In Weyl gravity, the 
metric is given by. 



-g{r)df + g{r)-^dr'^ + r'^^dO'^ + sin^ed(i?) 



9[r) 



1 

1 h 7r 

r 



kr 



(60) 
(61) 



Even though the context in which the Weyl gravity and the quintessence black hole 
are formulated are different, one can see that for k = and 7 = — c, the Weyl 
gravity geometry and the quintessence black hole geometry is the same!. Therefore, 
the bending angle obtained by Amore et.al[28] for the Weyl gravity can be modified 
to obtain the angle for the quintessence black hole as. 
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Note that here, To is the one given in eq.(56) for the quintessence black hole. 

When c — 0, one obtain the well known bending angle for the Schwarzschild black 
hole, as. 
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Here, corresponds to the one given in eq.(59). We will substitute the respective 
expressions for tq in terms of the impact paramter D for both black holes and plot 
the deflection angle as a function of the impact parameter D. 
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Figure 19. The graph shows the bending angle a for the quintessence black hole(dark 
curve) and the Schwarzschild black hole(dashed curve) as a function of the impact 
parameter D. Here, M = 1 and c = 0.05. = 6.8956 for the quintessence black 
hole and Dc = 5.19615 for the Schwarzschild black hole. 

From Fig. 19, the photons with the same impact parameter bends more around 
the quintessence black hole for small D and less for large D. We like to mention 
here that Liu et.al [31] did studied the gravitational frequency shift and deflection of 
light for the black holes surrounded by the quintessence for various values of the state 
parameter Uq. There, the function u was computed as a first order approximation 
where as here the exact solution is used. 



8 Conclusion 

We have studied the null geodesies of a black hole surrounded by dark energy. Here, 
we have chosen quintessence as the candidate for dark energy and the black hole 
solution studied was derived by Kiselev [14]. The equations for the geodesies were 
solved exactly for various values of energy and angular momentum of the photons. 
All possible motions are presented. The circular orbits are studied in detail and are 
shown to be unstable. We have also computed the Lyapunov exponent A which gives 
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the instability time scale for the unstable geodesies. It is shown that A is smaller for 
the dark energy black hole in comparison with the A for the unstable circular orbits 
for the Schwarzschild black hole. 

As an application of the photon motion studied here, we have studied the light de- 
flection for this particular black hole and have done comparisons with the Schwarzschild 
black hole. The photons with the same impact parameter bends more around the 
quintessence black hole for small D and less for large D. 

As an extension of this work, one could study the motion of particles around the 
charged version of the quintessence black hole. 
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